Low Temperature Specific Heat of Optimally Doped BaFe 2 _ a; TM :c As2 (TM = Co and 
Ni) Single Crystals: Constraint on the Pairing Gap 
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Low temperature specific heat has been measured in optimally doped and highly overdoped non- 
superconducting BaFe2-:rT]VLAs2 (TM = Co and Ni) single crystals. By using the data of the 
overdoped samples, we successfully removed the phonon contribution of the optimally doped ones, 
and derived the electronic specific heat coefficient y e . Remarkably, we found a continuing temper- 
ature dependent 7 e (T) which follows the quadratic relation 7 e = 70 + aT 2 in the low temperature 
limit. Together with the very small residual term 70, linear magnetic field dependence of -y e , it is 
concluded that there are either small segments of nodal lines, or point-like nodes in these samples. 

PACS numbers: 74.20.Rp, 74.70.Dd, 74.62.Dh, 65.40.Ba 
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The discovery of superconductivity above 50 K in iron 
pnictides has added a new member in the family of uncon- 
ventional high temperature superconductors. [l[ One of 
the key issues here is about the superconducting pairing 
mechanism. Theoretically it was suggested that the pair- 
ing may be established via inter-pocket scattering of elec- 
trons between the hole pockets (around T point) and elec- 
tron pockets (around M point), leading to the so-called 
S 1 * 1 pairing manner. 043 Experimental results about the 
pairing symmetry remain highly controversial leaving 
the perspectives ranging from S ++ -wave, to S and to 
d-wave.|6l4l5j Although evidence for a nodal gap has 
been accumulated in LaFePO@,0| and Ba(FeAsi_ x P x )2 
systems [13l|. in the charge doped 122, the experimental 
data point to the existence of isotropic gap s, especially 
in the optimally doped samples. [j[ ITlL Tl2j| The pene- 
tration depth measurements indicate a quadratic tem- 
perature dependence AA oc 1 — (T/T c ) 2 when the mea- 
suring inductive current is flowing along Fe As-plane. [T2]] 
And it becomes more linear like when the detecting cur- 
rent has a c-axis component, indicating the possibility 
of nodes along c-axis but nodeless along Fe As-plane .[la] 
This is further strengthened by the recent thermal con- 
ductivity measurements which suggest the possible ex- 
istence of a horizontal nodal line. 17] In this Letter, we 
present the data of low temperature specific heat (SH) in 
BaFe 2 -xTM a: As2 (TM = Co and Ni) single crystals. A 
unique quadratic temperature dependence 7 e = 7o + otT 2 
was found in the low-T limit. This is different from 
the expectation for a d-wave superconductor in the clean 
limit C e oc T 2 . Combining with the very small 70 and the 
linear field dependence of j e , it is tempting to conclude 
the existence of small segments of nodal lines or point- 
like nodes, being consistent with the recent five band 
tight binding calculations. (l8| 

The Co- and Ni-doped BaFeoAso single crystals were 
grown by the self-flux method[l9j,|20j]. In order to get the 
phonon contribution of SH, we also measured an over- 
doped non-superconducting sample as the reference for 



each kind of dopant. The typical dimension of the sin- 
gle crystals for specific heat measurements was 3x 3x 0.4 
mm 3 . The SH measurements were done with the thermal 
relaxation method on the Quantum Design instrument 
physical property measurement system (PPMS) with the 
temperature down to 1.8 K and magnetic field up to 9 
T. To improve the resolution, the field dependence of the 
sensor on the measuring chip was calibrated in advance. 

In Fig. 1(a), we show the temperature dependence of 
resistivity for four samples: BaFe2- :E Co a ;As2 with x — 
0.16 and 0.60, and BaFe 2 _ a ;Ni :E As 2 with x — 0.10 and 
0.30, respectively. The onset transition temperatures 
determined at p = 95%p„ for BaFei.84Coo.i6As2 and 
BaFe1.gNio.1As2 are 24.2 K and 20.5 K, respectively, 
which are close to the optimal doping points. The sharp 
transitions near T c in the dc magnetization and resistivity 
measurements indicate the good quality of the samples. 
No superconductivity was detected in BaFe1.4Coo.6As2 
and BaFe1.7Nio.3As2 down to 1.8 K. 

In Fig. 2, we present the SH of the superconducting 
and non-superconducting samples under zero field for 
each dopant. Sharp and clear SH anomalies near T c for 
the samples BaFei.84Coo.i6As2 and BaFe1.9Nio.1As2 can 
be seen. For all samples investigated here, the Schottky 
anomalies are negligible. For each dopant, one can see 
that the SH in the normal state for the superconduct- 
ing and non-superconducting samples share a very sim- 
ilar temperature dependence. This allows us to use the 
phonon part of SH of the non-superconducting sample 
as the reference, and subtract it safely for the supercon- 
ducting one. 

In order to extract the electronic SH of the supercon- 
ducting samples, people usually apply a magnetic field 
to destroy the superconductivity, and measure the SH 
of the normal state. For the iron-based superconduc- 
tors, however, the upper critical field is too high to be 
reached. I21B S o we employ the method suggested in the 
literatures |22l [23j . Here we take the overdoped non- 
superconducting sample as the references. Suppose that 
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FIG. 1: (color online) (a) Temperature dependence of re- 
sistivity of the optimally doped and highly overdoped non- 
superconducting BaFe 2 - : rTM E As2 (TM = Co and Ni) sin- 
gle crystals at zero field, (b) Temperature dependence of dc 
magnetization for the two superconducting samples. The en- 
larged view of resistivity for the two superconducting samples 
is shown in (c). 



the phonon contributions of the non-superconducting 
sample is C^ h (T), which can be obtained by subtracting 
the linear electronic term from the total SH, and that 
of the superconducting sample is C^ h (T), we naturally 



have C% h {T) = a ■ C^ h (b ■ T). Here a and b are fitting 
parameters which should be close to unity. Using a least- 
squares fit of our data above T c , we determined a and b to 
be 1.000, 1.013 for the Co-doped sample, and 0.98, 0.99 
for the Ni-doped one, indicating that the phonon con- 
tributions of the superconducting sample and the over- 
doped non-superconducting one are indeed very close to 
each other. Although the antiferromagnetic (AF) spin 
fluctuation is very strong in the optimally doped sample 
and negligible in the overdoped one 24] , the SH concerns 
mainly the Q = scattering and is sensitive only to the 
quasiparticle density of states (DOS) at Ep. The entropy 
conservation (not shown here) obtained by this treatment 
supports above argument. The electronic SH of the su- 
perconducting sample can then be obtained through 
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FIG. 2: (color online) Raw data of the temperature depen- 
dence of specific heat for the BaFe 2 -xTM :E As2 (TM = Co 
and Ni) single crystals. The enlarged views of the same data 
are presented as insets in (a) and (b), respectively. 



where Cf ot (T), Cf (T) are the total and the electronic SH 
of the superconducting sample, respectively. 

The obtained electronic SH for Co- and Ni-doped sam- 
ples are shown in Fig. 3. Surprisingly, we found a con- 
tinuing temperature dependent electronic SH coefficient 
7 e (T), which follows the quadratic relation j e — 70 + aT 2 
in the low-T limit (see the insets of Fig. 3). The parame- 
ters of 70 were determined from the fitting: 1.53 mJ/mol 
K 2 and 1.49 mJ/mol K 2 , for the superconducting Co- and 
Ni-doped samples, respectively, which are rather small 
compared to the values reported in literatures, indicating 
the cleanness of our samples. In order to get a compre- 
hensive understanding, we used the BCS formula to fit 
our data 
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where i e = le - 7o , C = V^+^W^), N F the 
DOS at the Fermi energy. We use four different gap 
structures to fit the data: single isotropic s-wave gap 
(Siso), single anisotropic gap with a d-wave feature (Sd- 
A = Ao|cos2#|), mixture of two: Si SO + Sd and two 
isotropic gaps Si SO i+Si S0 2- In the latter two cases, the 7 e 
was calculated through a linear combination of the two 
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FIG. 3: (color online) The electronic specific heat of the opti- 
mally doped superconducting samples (a) BaFei.84Coo.i6As2 
and (b) BaFe1.9Nio.1As2. These data were obtained by re- 
moving the phonon contribution (see text). The insets show 
the enlarged view of the low temperature specific heat, a 
quadratic temperature dependence j e = 70 + otT 2 can be 
clearly seen. Four different models: single isotropic s-wave 
gap (Siso), single anisotropic gap with a d-wave feature (Sd), 
mixture of two: Si so + Sd and Si so i + Si S0 2 were used to fit 
the data. 
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FIG. 4: (color online) Temperature dependence of the specific 
heat C/T in the low temperature region for BaFe1.84Coo.i6 As2 
and BaFe1.9Nio.1As2 with the magnetic field aligned with 
FeAs-plane and c-axis. 



TABLE I: Fitting parameters with different models for the 
Co-doped sample. 

model Ai(meV) fraction-1 A2(meV) fraction-2 



Si 

s d 

Siso ~\~ Sd 
Sisol H~ Siso2 



4.2 
5.4 
2.4 
1.0 



100% 
100% 
17% 
25% 



5.9 
4.25 



83% 
75% 



TABLE II: Fitting parameters with different models for the 
Ni-doped sample. 

model Ai(meV) fraction-1 A2(meV) fraction-2 



Siso 

Sd 

Siso ~\~ Sd 
Sisol 4~ Siso2 



3.1 
4.1 
2.0 
1.15 



100% 
100% 
29% 
29% 



4.7 
3.3 



71% 
71% 



components. The fitting curves are shown by the col- 
ored lines in Fig. 3, and the fitting parameters are listed 
in Table I and Table II, respectively. One can see that 
the single isotropic s-wave model can not fit the data at 
all, while the fittings with Sd, Si SO + Sd and Si so i + Si S0 2 
can roughly describe the data. However, if we explore the 
data in more detail, we can see that the fit of Si SO \ + Si S0 2 
model has an overall deviation to the data, and single Sd 
can roughly fit the data except in low-T region, where the 
fitting result is higher than the experimental data. The 
Siso + Sd model seems to give the best fitting in the whole 
temperature region. The multigap scenario (Si SO -\- Sd.^j 
with a very small s-wave gap of about 2 meV, which ac- 
counts for less than 30% of the total quasiparticle DOS, 
seems possible. The global fitting with Si SO +Sd, together 
with the relationship 7 e = 70 + aT 2 in the low-T limit, 
suggest that anisotropic gaps, possibly with nodes, may 



exist in these Fe As-based superconductors. 

To explore whether the quadratic temperature depen- 
dence 7 e = 70 + aT 2 is due to the presence of line nodes, 
we measured also the magnetic field dependence of j e . 
The data are plotted as C/T vs. T 2 in Fig. 4. One can see 
the roughly linear behavior in the low-T region. It is clear 
that the field induced enhancement of the electronic SH 
is larger when the field is along c-axis than along FeAs- 
plane, indicating that the upper critical field is higher 
along FeAs-plane than along c-axis. [25| By doing a linear 
extrapolation of the data in Fig. 4 to zero K, we get the 
field dependence of the electron SH coefficient j e (H), as 
shown in the upper part of Fig. 5. It is clear that 7e(-ff) 
increase linearly with magnetic field for the two samples 
in both alignments. This excludes both the vertical (like 
a d-wave) or horizontal line nodes (bottom-left picture in 
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FIG. 5: (color online) Magnetic field dependence of the elec- 
tronic specific heat for the samples BaFei.84Coo.i6As2 and 
BaFe1.9Nio.1As2 when the field is aligned along FeAs-plane 
and c-axis. A rough linear feature was observed in all four 
cases. The bottom part illustrates the possible gap nodes on 
the F Fermi surfaces: (left) horizontal nodal lines; (middle) 
small segments of nodal lines (shown by the red lines); (right) 
point nodes on the 3D Fermi pockets. 



Fig.5), since otherwise, as found in cupratc superconduc- 
tors, a square root relation "f e {H) oc v / 17[26] in the clean 
limit, or a curved feature "j e (H) oc H / H C 2log{BH c2 / H) 
with impurity scattering [27j should be observed. 

Although the quadratic temperature dependence j e — 
70 + aT 2 can also be attributed to the DOS induced by 
the impurity scattering within the S ± model 28 -3(|, we 
would argue that this is unlikely for the following rea- 
sons. Firstly, the exponent "2" here may be achieved 
in the case of strong pair breaking[3(|, which is in con- 
trast to the small value of 70 ~ 1.5 mJ/molK 2 (< 7.5% 
7 n ). Furthermore, according to Onari et al.,[29(, even 
with such a small amount impurities (assuming to induce 
the large momentum scattering), the superconductivity 
may not survive within the model. Secondly, it is 
not straightforward to understand the linear field depen- 
dence of 7 e (-ff) — 70 oc H with this picture. In order to 
reconcile the observations of a small 70, 7 e = 70 + olT 2 
and je(H) — 70 oc H, we propose that nodes, in the 
form of point-like or small segments, may exist in our 
present samples. As indicated by the five-band tight- 



binding calculations, [18J small segments of nodal lines (or 
called as accidental nodes) may exist on the hole-pockets 
near k z = ±n, which is depicted by the middle-bottom 
cartoon picture of Fig.5. An alternative interpretation 
would assume the extreme case of point-like nodes with 
a closed 3D Fermi pocket. It is natural that the pairing 
interaction by mediating the AF spin fluctuations should 
be very weak on the right top or bottom of these 3D 
pockets where the Fermi velocity is along c-axis (see the 
bottom- right picture of Fig.5), therefore the gap ampli- 
tude becomes negligible. Actually these two pictures can 
reconcile all the observations in this work. In a system 
with point nodes (or small segments of line nodes), it is 
anticipated that Np oc E 2 , since the Doppler shift energy 
of Eh oc VW, we have Np oc H. Meanwhile, the relation 
7e = 7o + aT 2 is also anticipated by the point-like nodes. 

In summary, the electronic specific heat on the opti- 
mally Ni- and Co-doped superconducting samples was 
derived by successfully removing the phonon contribu- 
tions. A quadratic temperature dependence of 7 e = 
70 + aT 2 in the low-T limit was discovered. The global 
temperature dependence of 7 e can be fitted with a two- 
component model with possible nodes. However the lin- 
ear field dependence 7 e (-ff) oc H observed with the mag- 
netic field along FeAs-plane and c-axis excludes the ex- 
istence of either vertical or horizontal line nodes. The 
results put strong constraint on the pairing gap and can 
be reconciled by the model of small segments of line nodes 
or point-like nodes. 
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